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Abstract 

The B~ -> D+K~ 7T decay is observed in a data sample corresponding to 3.0 fb 1 
of pp collision data recorded by the LHCb experiment during 2011 and 2012. Its 
branching fraction is measured to be B(B~ —y D + K~n~) = (7.31 ± 0.19 ± 0.22 ± 
0.39) x 10 -5 where the uncertainties are statistical, systematic and from the branching 
fraction of the normalisation channel B~ —> , respectively. An amplitude 

analysis of the resonant structure of the B~ —> D + K~-k~ decay is used to measure the 
contributions from quasi-two-body B~ -> Dq(2400)°K~, B~ —>• D2(2460)°I\~, and 
B~ —> D*j(2760)°I\~ decays, as well as from nonresonant sources. The D*j(2 760)° 
resonance is determined to have spin 1. 
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1 Introduction 


Excited charmed mesons are of great theoretical and experimental interest as they allow 
detailed studies of QCD in an interesting energy regime. Good progress has been achieved 
in identifying and measuring the parameters of the orbitally excited states, notably from 
Dalitz plot (DP) analyses of three-body B decays. Relevant examples include the studies of 
B~ —> D + n~ii~ 111 2] and B° —>■ D°7r + 7i~ |3 decays, which provide information on excited 


neutral and charged charmed mesons (collectively referred to as D** states), respectively. 
First results on excited charm-strange mesons have also recently been obtained with the 
DP analysis technique |4]-[6). Studies of prompt charm resonance production in e + e - 
and pp collisions [7,8| have revealed a number of additional high mass states. Most 
of these higher mass states are not yet confirmed by independent analyses, and their 
spectroscopic identification is unclear. Analyses of resonances produced directly from e + e~ 
and pp collisions do not allow determination of the quantum numbers of the produced 
states, but can distinguish whether or not they have natural spin parity (ie. J p in the 
series 0 + , 1 _ , 2 + ,...). The current experimental knowledge of the neutral D** states is 
summarised in Table [l] (here and throughout the paper, natural units with h = c = f 
are used). The £q(2400)°, £i( 2420)°, £((2430)° and £>2(2460)° mesons are generally 
understood to be the four orbitally excited (fP) states. The experimental situation as 
well as the spectroscopic identification of the heavier states is less clear. 

The B~ —>• D + K~n~ decay can be used to study neutral D** states. The D + K~n~ 
final state is expected to exhibit resonant structure only in the D + n~ channel, and unlike 
the Cabibbo-favoured final state does not contain any pair of identical particles. 

This simplifies the analysis of the contributing excited charm states, since partial wave 
analysis can be used to help determine the resonances that contribute. 

One further motivation to study B~ —>■ D + K~n~ decays is related to the measurement 

of the angle 7 of the unitarity triangle defined as 7 = arg [—V d V* b j ( V d V*,')]. where V 

--' 


are elements of the Cabibbo-Kobayashi-Maskawa (CKM) quark mixing matrix 10,11 


Table 1: Measured properties of neutral D** states. Where more than one uncertainty is given, 
the first is statistical and the others systematic. 


Resonance 

Mass 

Width 

Ref. 


(MeV) 

(MeV) 



£>g(2400)° 

2318 ± 29 

267 ± 40 

0 + 

9 

£i(2420)° 

2421.4 ±0.6 

27.4 ±2.5 

1 + 

9 

£>((2430)° 

2427 ± 26 ± 20 ± 15 

384 +*° 7 ± 24 ± 70 

1 + 

T 

£>2(2460)° 

2462.6 ± 0.6 

49.0 ± 1.3 

2 + 

9 

£>*(2600) 

2608.7 ±2.4 ±2.5 

93 ± 6 ± 13 

natural 

7 

£>*(2650) 

2649.2 ±3.5 ±3.5 

140 ± 17 ± 19 

natural 

8 

£*(2760) 

2763.3 ±2.3 ±2.3 

60.9 ±5.1 ±3.6 

natural 

7 

£*(2760) 

2760.1 ±1.1 ±3.7 

74.4 ± 3.4 ± 19.1 

natural 
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One of the most powerful methods to determine 7 uses B 
neutral D meson decaying to CP eigenstates 


12,13 


DK decays, with the 
The sensitivity to 7 arises due to 


the interference of amplitudes proportional to the CKM matrix elements V ub and V c 


cbi 


associated with D° and D° production respectively. However, a challenge for such methods 
is to determine the ratio of magnitudes of the two amplitudes, Tb, that must be known 
to extract 7 . This is usually handled by including D meson decays to additional final 
states in the analysis. By contrast, in B~ —» D**K~ decays the efficiency-corrected ratio 
of yields of B~ — » D**K~ — > D~tt + K~ and B~ — » D**K~ —> D + tt~K^ decays gives r B 
directly 14 . The decay B~ —y D**K~ —y Dti 0 K~ where the D meson is reconstructed 
in CP eigenstates can be used to search for CP violation driven by 7 . Measurement of 
the first two of these processes would therefore provide knowledge of r B in B~ —y D**K~ 
decays, indicating whether or not a competitive measurement of 7 can be made with this 
approach. 

In this paper, the B~ —> D + K~7r~ decay is studied for the first time, with the 
D + meson reconstructed through the K~n + n + decay mode. The inclusion of charge 
conjugate processes is implied. The topologically similar B~ —>• D + 'ji~'k~ decay is used as 
a control channel and for normalisation of the branching fraction measurement. A large 
B~ — >■ D + K~tt~ signal yield is found, corresponding to a clear first observation of the 
decay, and allowing investigation of the DP structure of the decay. The amplitude analysis 
allows studies of known resonances, searches for higher mass states and measurement of 
the properties, including the quantum numbers, of any resonances that are observed. The 
analysis is based on a data sample corresponding to an integrated luminosity of 3.0 ftP 1 of 
pp collision data collected with the LHCb detector, approximately one third of which was 
collected during 2011 when the collision centre-of-mass energy was a/s = 7TeV and the 
rest during 2012 with a/s = 8 TeV. 

The paper is organised as follows. A brief description of the LHCb detector as well as 
reconstruction and simulation software is given in Sec. [2] The selection of signal candidates 
is described in Sec. [3| and the branching fraction measurement is presented in Sec. |4j 
Studies of the backgrounds and the fit to the B candidate invariant mass distribution 
are in Sec. HD with studies of the signal efficiency and a definition of the square Dalitz 
plot (SDP) in Sec. 4.2 Systematic uncertainties on, and the results for, the branching 


fraction are discussed in Secs. EH and |4.4| respectively. A study of the angular moments 
of B~ —> D + K~tt~ decays is given in Sec. [5j with results used to guide the Dalitz plot- 
analysis that follows. An overview of the Dalitz plot analysis formalism is given in Sec. [6], 
and details of the implementation of the amplitude analysis are presented in Sec. [7J The 
evaluation of systematic uncertainties is described in Sec. [8j The results and a summary 
are given in Sec. [9j 


2 LHCb detector 

The LHCb detector [l5p6] is a single-arm forward spectrometer covering the pseudorapidity 
range 2 < rj < 5, designed for the study of particles containing b or c quarks. The detector 
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includes a high-precision tracking system consisting of a silicon-strip vertex detector 0 
surrounding the pp interaction region, a large-area silicon-strip detector located upstream 
of a dipole magnet with a bending power of about 4 Tm, and three stations of silicon-strip 
detectors and straw drift tubes 18 placed downstream of the magnet. The polarity of 


the dipole magnet is reversed periodically throughout data-taking. The tracking system 
provides a measurement of momentum, p, of charged particles with a relative uncertainty 
that varies from 0.5% at low momentum to 1.0% at 200 GeV. The minimum distance of a 
track to a primary vertex, the impact parameter (IP), is measured with a resolution of 
(15 + 29/pr) pm, where Pt is the component of the momentum transverse to the beam, 
in GeV. Different types of charged hadrons are distinguished using information from 
two ring-imaging Cherenkov detectors 19 . Photon, electron and hadron candidates are 


identified by a calorimeter system consisting of scintillating-pad and preshower detectors, 
an electromagnetic calorimeter and a hadronic calorimeter. Muons are identified by a 
system composed of alternating layers of iron and multiwire proportional chambers 


20 


The trigger 21 consists of a hardware stage, based on information from the calorimeter 
and muon systems, followed by a software stage, in which all tracks with p t > 500 (300) MeV 
are reconstructed for data collected in 2011 (2012). The software trigger line used in 
the analysis reported in this paper requires a two-, three- or four-track secondary vertex 
with significant displacement from the primary pp interaction vertices (PVs). At least 
one charged particle must have p\ > 1.7 GeV and be inconsistent with originating from 
the PV. A multivariate algorithm [221 is used for the identification of secondary vertices 
consistent with the decay of a b hadron. 

In the offline selection, the objects that fired the trigger are associated with recon¬ 
structed particles. Selection requirements can therefore be made not only on the trigger 
line that fired, but on whether the decision was due to the signal candidate, other particles 
produced in the pp collision, or a combination of both. Signal candidates are accepted 
offline if one of the final state particles created a cluster in the hadronic calorimeter with 
sufficient transverse energy to fire the hardware trigger. These candidates are referred 
to as “triggered on signal” or TOS. Events that are triggered at the hardware level by 
another particle in the event, referred to as “triggered independent of signal” or TIS, are 
also retained. After all selection requirements are imposed, 57 % of events in the sample 
were triggered by the decay products of the signal candidate (TOS), while the remainder 
were triggered only by another particle in the event (TIS-only). 

Simulated events are used to characterise the detector response to signal and certain 


types of background events. In the simulation, pp collisions are generated using Pythia 23 
with a specific LHCb conhguration 24 . Decays of hadronic particles are described 

The 


in which final state radiation is generated using PHOTOS 26 


by EvtGen 

interaction of the generated particles with the detector and its response are implemented 
using the Geant 4 toolkit 


27 as described in Ref. 28 
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3 Selection requirements 


Most selection requirements are optimised using the B~ —y control channel. 

Loose initial selection requirements on the quality of the tracks combined to form the 
B candidate, as well as on their p. p T and y 2 P , are applied to obtain a visible peak in 
the invariant mass distribution. The Xip is the difference between the y 2 of the PV 
reconstruction with and without the considered particle. Only candidates with invariant 
mass in the range 1770 < m(K~'K + n + ) < 1968 MeV are retained. Further requirements 
are imposed on the vertex quality (y 2 tx ) an d flight distance from the associated PV of 
the B and D candidates. The B candidate must also satisfy requirements on its invariant 
mass and on the cosine of the angle between the momentum vector and the line joining 
the PV under consideration to the B vertex (cos#di r ). The initial selection requirements 
are found to be about 90 % efficient on simulated signal decays. 

Two neural networks 29 are used to further separate signal from background. The 
first is designed to separate candidates that contain real D + —> K ~ti + 'k + decays from 
those that do not; the second separates B~ —> D + ti~7t~ signal decays from background 
combinations. Both networks are trained using the D + n~'Ji~ control channel, where the 
sPlot technique 30 is used to statistically separate B~ —$■ signal decays from 


background combinations using the D ( B ) candidate mass as the discriminating variable 
for the first (second) network. The first network takes as input properties of the D 
candidate and its daughter tracks, including information about kinematics, track and 
vertex quality. The second uses a total of 27 input variables. They include the y 2 P °f the 
two “bachelor” pions (i.e. pions that originate directly from the B decay) and properties 
of the D candidate including its y 2 P , y 2 tx , cos#^, the output of the D neural network 
and the square of the flight distance divided by its uncertainty squared (Xfl ight ). Variables 
associated to the B candidate are also used, including p T , y 2 p , y 2 tx , Xfl ight and cos#^. 
The asymmetry and track multiplicity in a cone with half-angle of 1.5 units of the plane 
of pseudorapidity and azimuthal angle (measured in radians) around the B candidate 
flight direction 31 , which contain information about the isolation of the B candidate from 


the rest of the event, are also used in the network. The neural network input quantities 
depend only weakly on the kinematics of the B decay. A requirement is imposed on the 
second neural network output that reduces the combinatorial background by an order of 
magnitude while retaining about 75 % of the signal. 

The selection criteria for the B~ —> D + K~n~ and B~ —> D + candidates are 
identical except for the particle identification (PID) requirement on the bachelor track 
that differs between the two modes. All five final state particles for each decay mode 
have PID criteria applied to preferentially select either pions or kaons. Tight requirements 
are placed on the higher momentum pion from the D + decay and on the bachelor kaon 
in B~ —>- D + K~n~ to suppress backgrounds from Df K~K + ir + and B~ —>- D +/ k~ti^ 
decays, respectively. The combined efficiency of the PID requirements on the five final state 
tracks is around 70% for B~ —>• D + n~7r^ decays and around 40% for B~ -P D + K~tt~ 
decays. The PID efficiency depends on the kinematics of the tracks, as described in detail 
in Sec. 4.2, and is determined using samples of D° -P K~tt + decays selected in data by 
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exploiting the kinematics of the D* + —> D 0 tc + decay chain to obtain clean samples without 
using the PID information. 


To improve the B candidate invariant mass resolution, track momenta are scaled 32 33 


with calibration parameters determined by matching the measured peak of the J/^> —> /i + /i~ 
decay to the known J/ip mass |9|. Furthermore, a fit to the kinematics and topology of the 
decay chain 34 is used to adjust the four-momenta of the tracks from the D candidate so 


that their combined invariant mass matches the world average value for the D + meson |9|. 
An additional B mass constraint is applied in the calculation of the variables that are 
used in the Dalitz plot fit. 

To remove potential background from misreconstructed A f t decays, candidates are 
rejected if the invariant mass of the D candidate lies in the range 2280-2300 MeV when 
the proton mass hypothesis is applied to the low momentum pion track. Possible back¬ 
grounds from B~ meson decays without an intermediate charm meson are suppressed by 
the requirement on the output value from the first neural network, and any surviving 
background of this type is removed by requiring that the D candidate vertex is displaced 
by at least 1 mm from the B decay vertex. The efficiency of this requirement is about 
85%. 


Signal candidates are retained for further analysis if they have an invariant mass in 
the range 5100-5800 MeV. After all selection requirements are applied, fewer than 1% of 
events with one candidate also contain a second candidate. Such multiple candidates are 
retained and treated in the same manner as other candidates; the associated systematic 
uncertainty is negligible. 


4 Branching fraction determination 


The ratio of branching fractions is calculated from the signal yields with event-by-event 
efficiency corrections applied as a function of square Dalitz plot position. The calculation 

is 


B{B~ -> D+K~ 7T“) _ N COII (B- ->■ D+K~ir~) 
B{B~ D+ 7T-7T-) “ N C ™(B- ->■ D+7T-7T-) ’ 


where iV corr = W % /is the efficiency-corrected yield. The index i sums over all 

candidates in the data sample and W t is the signal weight for each candidate, which is 
determined from the fits described in Sec. 4. 1| and shown in Figs. [T] and [2j using the 
sPlot technique 30 . Each fit is performed simultaneously to decays in the TOS and 


TIS-only categories. The efficiency of candidate i , e*, is obtained separately for each trigger 
subsample as described in Sec. 4.2| 


4.1 Determination of signal and background yields 

The candidates that survive the selection requirements are comprised of signal decays 
and various categories of background. Combinatorial background arises from random 
combinations of tracks (possibly including a real D + —y K~tt + tt + decay). Partially 
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reconstructed backgrounds originate from b hadron decays with additional particles that 
are not part of the reconstructed decay chain. Misidentihed decays also originate from 
b hadron decays, but where one of the final state particles has been incorrectly identified 
( e.g. a pion as a kaon). The signal (normalisation channel) and background yields are 
obtained from unbinned maximum likelihood fits to the D + K~n~ (D + 'n~n~) invariant 
mass distributions. 

Both the B~ —s- D + K~n~ and B~ —> D + 7i~7i~ signal shapes are modelled by the 
sum of two Crystal Ball (CB) functions 35 with a common mean and tails on opposite 
sides, where the high-mass tail accounts for non-Gaussian reconstruction effects. The 
ratio of widths of the CB shapes and the relative normalisation of the narrower CB shape 
are constrained within their uncertainties to the values found in fits to simulated signal 
samples. The tail parameters of the CB shapes are also fixed to those found in simulation. 

The combinatorial backgrounds in both D + K~n~ and D +/ n~7i~ samples are modelled 
with linear functions; the slope of this function is allowed to differ between the two trigger 
subsamples. The decay B~ —> D* + K~n~ is a partially reconstructed background for 
D + K~tt~ candidates, where the D* + decays to either D + 7 or D + 7r° and the neutral particle 
is not reconstructed. Similarly the decay B~ —s- D* + 7r~7i~ forms a partially reconstructed 
background to the D + h~ti~ final state. These are modelled with non-parametric shapes 
determined from simulated samples. The shapes are characterised by a sharp edge around 
100 MeV below the B peak, where the exact position of the edge depends on properties 
of the decay including the D* + polarisation. The fit quality improves when the shape is 
allowed to be offset by a small shift that is determined from the data. 

Most potential sources of misidentihed backgrounds have broad B candidate invariant 
mass distributions, and hence are absorbed in the combinatorial background component 
in the fit. The decays B —> Th*) + 7 r _ 7 r~ and B~ —> DfK~TT~, however, give distinctive 
shapes in the mass distribution of D + K~n~ candidates. For D +/ k~ti~ candidates the only 
significant misidentihed background contribution is from B~ —>- D^ + K~ decays. The 
misidentihed background shapes are also modelled with non-parametric shapes determined 
from simulated samples. 

The simulated samples used to obtain signal and background shapes are generated 
with hat distributions in the phase space of their SDPs. For B~ —>■ D + 'k~h~ and B~ —y 
D* + n~TT~ decays, accurate models of the distributions across the SDP are known jTJ,|2j , so 
the simulated samples are reweighted using the B~ —>■ ~ data sample; this affects 

the shape of the misidentihed background component in the ht to the D + K~ti + sample. 
Additionally, the D + and D* + portions of this background are combined according to 
their known branching fractions. All of the shapes, except for that of the combinatorial 
background, are common between the two trigger subsamples in each ht, but the signal and 
background yields in the subsamples are independent. In total there are 15 free parameters 
in the ht to the D + h~ti~ sample: yields in each subsample for signal, combinatorial, 
B~ D&+K- "7T and B —> D* + 7 r 7 t backgrounds; the combinatorial slope in each 
subsample; the double CB peak position, the width of the narrower CB, the ratio of CB 
widths and the fraction of entries in the narrower CB shape; and the shift parameter of 
the partially reconstructed background. The result of the ht is shown in Fig. [l] 
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Figure 1: Results of the fit to the B~ —> candidate invariant mass distribution for the 

(left) TOS and (right) TIS-only subsamples. Data points are shown in black, the full fitted model 
as solid blue lines and the components as shown in the legend. 


Table 2: Yields of the various components in the fit to B —> D + n n candidate invariant mass 
distribution. 


Component 

TOS 

TIS-only 

N(B~ —> D+7T-7T-) 

29 190 ± 204 

19 416 ± 159 

N(B~ -)■ DW+K-tt-) 

807 ± 123 

401 ± 84 

N(B~ D* + n~7r~) 

12120 ± 115 

8551 ± 96 

fV(comb. bkgd.) 

784 ± 54 

746 ± 47 


for both trigger subsamples and gives a combined signal yield of approximately 49 000 
decays. Component yields are given in Table [2] 

There are a total of 17 free parameters in the fit to the D + K~tt~ sample: yields 
in each subsample for signal, combinatorial, B~ —>■ D* + K~n~, B~ —>■ and 

B~ DW+ 7r 7T backgrounds; the combinatorial slope in each subsample; the same 

signal shape parameters as for the D +/ k~ti~ fit; and the shift parameter of the partially 
reconstructed background. Figure [2] shows the result of the D + K~tt~ fit for the two trigger 
subsamples that yield a total of approximately 2000 B~ —>• D + K~tt~ decays. The yields 
of all fit components are shown in Table [3| The statistical signal significance, estimated 
in the conventional way from the change in negative log-likelihood from the fit when the 
signal component is removed, is in excess of 60 standard deviations (a). 

4.2 Signal efficiency 

Since both B~ —>- D + K~ty~ and B —> D + 7i~n~ decays have non-trivial DP distributions, 
it is necessary to understand the variation of the efficiency across the phase space. Since, 
moreover, the efficiency variation tends to be strongest close to the kinematic boundaries 
of the conventional Dalitz plot, it is convenient to model these effects in terms of the SDP 
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Figure 2: Results of the fit to the B —> D + K ~7T candidate invariant mass distribution for 
the (left) TOS and (right) TIS-only subsamples. Data points are shown in black, the full fitted 
model as solid blue lines and the components as shown in the legend. 

Table 3: Yields of the various components in the fit to B~ —> D + K ~candidate invariant 
mass distribution. 


Component 
N{B~ ->• D + K~tt-) 
N(B~ -> £>(*)+7T-7T-) 
N(B~ -> D+K~n-) 
N(B~ —> D*+K~ 7T-) 
IV (comb. bkgd.) 


TOS TIS-only 
1112 ±37 891 ±32 
114 ±34 23 ±27 

69 ± 17 40 ± 15 

518 ±26 361 ±21 

238 ± 38 253 ± 36 


defined by variables m' and 9' which are valid in the range 0 to 1 and are given for the 
D + K~7i~ case by 


m 


1 

— arccos 
7r 


m(D + 7T ) — 

rry-. max _ rr) min 

lib r~»_ I lib y—\_1_ 

7r D +7r 



and 9' = -9(D + ir~), (2) 

7T 


where = mg- - and m 1 ^ 1 + 7r _ = mn+ ± m n - are the kinematic boundaries 

of m(D + 7i~) allowed in the B~ —>• D + K~n~ decay and 9{D + n~) is the hclicity angle of 
the D + 7r~ system (the angle between the K~ and the D + meson momenta in the D + n~ 
rest frame). For the ~ case, m' and 9' are defined in terms of the 7r - 7r - mass 

and hclicity angle, respectively, since with this choice only the region of the SDP with 
9'(tt~71~) < 0.5 is populated due to the symmetry of the two pions in the final state. 

Efficiency variation across the SDP is caused by the detector acceptance and by trigger, 
selection and PID requirements. The efficiency variation is evaluated for both D + K~ir~ 
and final states with simulated samples generated uniformly over the SDP. Data- 

driven corrections are applied to correct for known differences between data and simulation 
in the tracking, trigger and PID efficiencies, using identical methods to those described in 
Ref. [5]. The efficiency functions are fitted with two-dimensional cubic splines to smooth 
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Figure 3: Signal efficiency across the SDP for (left) TOS and (right) TIS-only B —> D + K n 
decays. The relative uncertainty at each point is typically 5 %. 


out statistical fluctuations due to limited sample size. 

The efficiency is studied separately for the TOS and TIS-only categories. The efficiency 
maps for each trigger subsample are shown for B~ —> D + K ~decays in Fig. [ 3 ] Regions 
of relatively high efficiency are seen where all decay products have comparable momentum 
in the B rest frame; the efficiency drops sharply in regions with a low momentum bachelor 
track due to geometrical effects. The efficiency maps are used to calculate the ratio of 
branching fractions and also as inputs to the D + K~n~ Dalitz plot fit. 

4.3 Systematic uncertainties 

Table [4] summarises the systematic uncertainties on the measurement of the ratio of 
branching fractions. Selection effects cancel in the ratio of branching fractions, except for 
inefficiency due to the T+ veto. The invariant mass fits are repeated both with a wider 
veto (2270-2310 MeV) and with no veto, and changes in the yields are used to assign a 
relative systematic uncertainty of 0.2%. 

To estimate the uncertainty arising from the choice of invariant mass fit- model, 
the D + K~n~ mass fit is varied by replacing the signal shape with the sum of two 
bifurcated Gaussian functions, removing the smoothing of the non-parametric functions, 
using exponential and second-order polynomial functions to describe the combinatorial 
background, varying fixed parameters within their uncertainties and varying the binning of 
histograms used to reweight the simulated background samples. For the D + n~7r - fit the 
same variations are made. The relative changes in the yields are summed in quadrature to 
give a relative systematic uncertainty on the ratio of branching fractions of 2.0%. 

The systematic uncertainty due to PID is estimated by accounting for three sources: 
the intrinsic uncertainty of the calibration (1.0%); possible differences in the kinematics 
of tracks in simulated samples, used to reweight the calibration data samples, to those 
in the data (1.7%); the granularity of the binning in the reweighting procedure (0.7%). 
Combining these in quadrature, the total relative systematic uncertainty from PID is 
2 . 1 %. 
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Table 4: Relative systematic uncertainties on the measurement of the ratio of branching fractions 
for B~ —> D + K~r~ and B~ —> D + r~r~ decays. 


Source 

Uncertainty (%) 

T+ veto 

0.2 

Fit model 

2.0 

Particle identification 

2.1 

Efficiency modelling 

0.8 

Total 

3.0 


The bins of the efficiency maps are varied within uncertainties to make fOO new 
efficiency maps, for both D + K~tt~ and D + r~r~ modes. The efficiency-corrected yields 
are evaluated for each new map and their distributions are fitted with Gaussian functions. 
The widths of these are used to assign a relative systematic uncertainty on the ratio of 
branching fractions of 0.8%. 

A number of additional cross-checks are performed to test the branching fraction result. 
The neural network and PID requirements are both tightened and loosened. The data 
sample is divided by dipole magnet polarity and year of data taking. The branching 
fraction is also calculated separately for TOS and TIS-only events. All cross-checks give 
consistent results. 

4.4 Results 

The ratio of branching fractions is found to be 
B(B~ -)• D+K-ir-) 

-—-f = 0.0720 ± 0.0019 ± 0.0021, 

B(B~ —> D + R~R~) 

where the first uncertainty is statistical and the second systematic. The statistical 
uncertainty includes contributions from the event weighting used in Eq. (JTj) and from the 
shape parameters that are allowed to vary in the fit [361. The world average value of 
B(B~ — > D + r~r~) = (1.07T0.05) x 10~ 3 |9| assumes that B + B~ and B°B° are produced 
equally in the decay of the T(4S') resonance. Using T(T(4S') —>■ B + B~)/T(T(4S) — > 
B°B° ) = 1.055 ± 0.025 |9| gives a corrected value of B(B~ —> D + r~r~) = (1.01 ± 0.05) x 
10~ 3 . This allows the branching fraction of B~ —>■ D + K~r~ decays to be determined as 

B{B~ ->■ D + K~r~) = (7.31 ± 0.19 ± 0.22 ± 0.39) x 10“ 5 , 

where the third uncertainty is from B(B~ —y D + r~r~). This measurement represents the 
first observation of the B~ —> D + K~r~ decay. 

5 Study of angular moments 

To investigate which amplitudes should be included in the DP analysis of B~ —> D + K~r~ 
decays, a study of its angular moments is performed. Such an analysis is particularly 
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useful for B~ —> D + K~n~ decays because resonant contributions are only expected to 
appear in the D + n~ combination, and therefore the distributions should be free of effects 
from reflections that make them more difficult to interpret. 

The analysis is performed by calculating moments from the Legendre polynomials 
Pl of order up to 2 J max , where J max is the maximum spin of the resonances considered. 
Each candidate is weighted according to its value of P h (cos 6(D + n~)) with an efficiency 
correction applied, and background contributions subtracted. The results for J max = 3 are 
shown in Fig. [I] for the D + n~ invariant mass range 2.0-3.0 GeV. The distributions of (P$) 
and ( Pq ) are compatible with being flat, which implies that there are no significant spin 3 
contributions. Considering only contributions up to spin 2, the following expressions are 
used to interpret Fig. [4] 


(Po) OC l^o | 2 + |^11 2 + | ^21 2 ) 

2 4 

(Pi) cx —= \h Q \ \hi\ cos (h 0 - hi) H —-= \hi\ \h 2 \ cos (hi - h 2 ) 


V3 


VT5 


(P 2 ) cx |ho| \h 2 \ cos (ho — h 2 ) + - |hi | " + - \h 2 \ 
a/5 5 7 


(F 3 ) |hi| |h 2 1 cos (hi — h 2 ) , 

(P 4 ) oc| |h 2 | 2 , 


(3) 

(4) 

(5) 

( 6 ) 
(7) 


where S-, P- and D-wave contributions are denoted by amplitudes hje l5j (j = 0,1,2 
respectively). The Ph;(2460) 0 resonance is clearly seen in the (P 4 ) distribution of Fig.jdjh). 
The distribution of (P 3 ) shows interference between spin 1 and 2 contributions, indicating 
the presence of a broad, possibly nonresonant, spin 1 contribution at low m(D + 7T~). 
The difference in shape between (Pi) and (P 3 ) shows interference between spin 1 and 0 
indicating that a broad spin 0 component is similarly needed. 


6 Dalitz plot analysis formalism 


A Dalitz plot 1371 is a representation of the phase-space for a three-body decay in terms 
of two of the three possible two-body invariant mass squared combinations. In /F —y 
D + K~7 r~ decays, resonances are expected in the m 2 {D + Ti~) combination, therefore this 
and m 2 (D + K~) are chosen to define the DP axes. For a fixed B mass, all other relevant 
kinematic quantities can be calculated from these two invariant mass squared combinations. 

The complex decay amplitude is described using the isobar approach 38-40 , where 


the total amplitude is calculated as a coherent sum of amplitudes from resonant and 
nonresonant intermediate processes. The total amplitude is then given by 


N 


A(m 2 {D + 7 T ),m 2 (D + K )) = ^^CjFj (m 2 ( y D + 7T ),m 2 (D + K )) , 


( 8 ) 


3 = 1 
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Figure 4: The first seven Legendre-polynomial weighted moments for background-subtracted 
and efficiency-corrected B~ —> D + K~tt~ data (black points) as a function of m(D + n ~) in the 
range 2.0-3.0 GeV. Candidates from both TOS and TIS-only subsamples are included. The blue 
line shows the result of the DP fit described in Sec. [7l 
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where Cj are complex coefficients giving the relative contribution of each intermediate 
process. The Fj (m 2 , m 2 (D + K~)) terms contain the resonance dynamics, which 
are composed of several terms and are normalised such that the integral of the squared 
magnitude over the DP is unity for each term. For a D + n~ resonance 

F (m 2 (D + 7r-),m 2 (D + K-)) = R(m(D + 7r-)) x X(\p\r BW ) x X(\q\r BW ) xT(p,q), (9) 


where the functions R, X and T are described below, and p and q are the bachelor particle 
momentum and the momentum of one of the resonance daughters, respectively, both 
evaluated in the D + tt~ rest frame. 

The X(z) terms, where z = \q | rew or \p\ tbw, are Blatt-Weisskopf barrier factors 41 
with barrier radius r B w> and are given by 


L = 0 : X(z) 
L = 1 : X(z) 
L = 2 : X(z) 
L = 3 : X(z) 


l±fo 

1+z 2 


z o + 3z o 
z 4 +3z 2 +9 ’ 

z$+6z$+45z$+225 
z 6 +6z 4 +45z 2 +225 ’ 


( 10 ) 


where Zq is the value of z when the invariant mass is equal to the pole mass of the resonance 
and L is the spin of the resonance. For a D + n~ resonance, since the B~ meson has zero 
spin, L is also the orbital angular momentum between the_resonance and the kaon. The 
barrier radius, tew, is taken to be 4.0 GeV” 1 ~ 0.8 fin (5,42 for all resonances. 

The terms T(p,q) describe the angular probability distribution and are given in the 
Zenrach tensor formalism |43,|44| by 


L = 0 

T(p,q) = 

1, 

L = 1 

T{p,q) = 

— 2 p-q, 

L = 2 

T{p,q) = 

l[3(p-q) 2 -(\p\\q\) 2 } , 

L = 3 

T{p,q) = 

- it [5(p • q) 3 - 3 (p- q)(\p\\q\) 2 


( 11 ) 


which are proportional to the Legendre polynomials, Pl(x), where x is the cosine of the 
angle between p and q (referred to as the hclicity angle). 

The function R (m(D + n~)) of Eq. Q is the mass lineshape. The resonant contributions 
considered in the DP model are described by the relativistic Breit-Wigner (RBW) function 

= (mg-m 2 ) -im 0 r(m) ’ ^ 

where the mass-dependent decay width is 

/ \ 2L+1 

r(m) = r 0 ( — ) (— ) X 2 (qr BW ) , (13) 

\q 0 J \mJ 


where qo is the value of q — \q\ for m — m o. Virtual contributions, from resonances with 
pole masses outside the kinematically accessible region of the phase space, can also be 
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modelled by this shape with one modification: the pole mass m 0 is replaced with 
a mass in the kinematically allowed region, in the calculation of the parameter q 0 . This 
effective mass is defined by the ad hoc formula [5] 


mf(m 0 ) = rn mm + (m max - rn mm ) 1 + tanh 


m 0 - 


a +m n 


m 1 


— rrv 


(14) 


where m max and rn mm are the upper and lower limits of the kinematically allowed range, 
respectively. For virtual contributions, only the tail of the RBW function enters the Dalitz 
plot. 

Given the large available phase-space in the B decay, it is possible to have nonresonant 
amplitudes (be. contributions that are not from any known resonance, including virtual 
states) that vary across the Dalitz plot. A model that has been found to describe well 
nonresonant contributions in several B decay DP analyses is an exponential form factor 
(EFF) (4g, 

R(m) = e ~ am2 , (15) 

where m is a two-body (in this case Dn) invariant mass and a is a shape parameter that 
must be determined from the data. 

Neglecting reconstruction effects, the DP probability density function would be 


V phys (m 2 (D + 7T ),m 2 (D + K )) 


\A(m 2 (D + n-) 1 m 2 (D + K-)) \ 2 
//dp l^l 2 dm 2 (D + n~) dm 2 (D + K~) ’ 


(16) 


where the dependence of A on the DP position has been suppressed in the denominator 
for brevity. The complex coefficients, given by Cj in Eq. ([8]), are the primary results of 
most Dalitz plot analyses. However, these depend on the choice of normalisation, phase 
convention and amplitude formalism in each analysis. Fit fractions and interference fit- 
fractions are also reported as these provide a convention-independent method to allow 
meaningful comparisons of results. The fit fraction is defined as the integral of the 
amplitude for a single component squared divided by that of the coherent matrix element 
squared for the complete Dalitz plot, 


FFi 


//dp I c i F o (m 2 (D + 7r-),m 2 (D+K-))\ 2 dm 2 (D + n~) dm 2 (D + K~) 
//dp l*4| 2 dm 2 (D + 7T _ ) dm 2 (D + K~) 


(17) 


The fit fractions do not necessarily sum to unity due to the potential presence of net- 
constructive or destructive interference, described by interference fit fractions defined for 
i < j only by 


FF ■ ■ — 
rr y 


/ f Dp 2 77e [ Cl (f 3 F t F*] dm 2 (D + -K~) dm 2 {D+K- 
//dp 1-4/ dm 2 (D + 7Tdm 2 (D + K~) 


( 18 ) 


where the dependence of F} ; and A on the DP position has been omitted. 
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Figure 5: Distribution of B —> D + K n candidates in the signal region over (left) the DP and 
(right) the SDP. Candidates from both TOS and TIS-only subsamples are included. 


7 Dalitz plot fit 


The LAURA±± (46 package is used to perform the Dalitz plot fit, with the two trigger 


subsamples fitted simultaneously using the Jfit method 47 . The two subsamples have 


separate signal and background yields, efficiency maps and background SDP distributions, 
but all parameters of the signal model are common. The likelihood function that is used is 


N c 




Y NkV k (m;lm?(D + K-)) 


(19) 


where the index i runs over N c candidates, while k distinguishes the signal and background 
components with TV*, the yield in each component. The probability density function for 
signal events, 77 


Slg , is given by Eq. (16) where the \A(m' 2 (D + n ),m 2 (D + K )) | 2 terms 


are multiplied by the efficiency function described in Sec. El The mass resolution is 
approximately 2.4 MeV, which is much lower than the width of the narrowest contribution 
to the Dalitz plot (~ 50 MeV); therefore, this has negligible effect on the likelihood and is 
not considered further. 

The signal and background yields that enter the Dalitz plot fit are taken from the mass 
fit described in Sec. |4.1| Only candidates in the signal region, defined as ±2.5 a around the 
B signal peak, where a is the width of the peak, are used in the Dalitz plot fit. Within this 
region, in the TOS subsample the result of the B candidate invariant mass fit corresponds 
to yields of 1060 ± 35, 37 ± 6, 26 ± 8 and 16 ± 4 in the signal, combinatorial background, 
£>(*)+ 7 r 7T and D+K it components, respectively. The equivalent yields in the TIS-only 
subsample are 849 ± 30, 39 ± 6, 5 ± 5 and 9 ± 3 candidates. The contribution from 
D* + K~tt~ decays is negligible in the signal window. The distributions of the candidates 
in the signal region over the DP and SDP are shown in Fig. [5] 

The SDP distributions of the Z^*) + 7r _ 7r _ and DfK~7T~ background sources are ob¬ 
tained from simulated samples using the same procedures as described for their invariant 
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Figure 6: Square Dalitz plot distributions used in the Dalitz plot fit for (top) combinatorial back¬ 
ground, (middle) B~ —> D ^*> + decays and (bottom) B~ —> Df K~tt~ decays. Candidates 
the TOS (TIS-only) subsamples are shown in the left (right) column. 


mass distributions in Sec. 4.1 


The distribution of combinatorial background events is mod- 
candidates in the sideband high-mass range 5500-5800 MeV, 
in this region subtracted. The dependence of the SDP 


elled by considering D + K~tt 
with contributions from 
distribution on B candidate mass was investigated and found to be negligible. The SDP 
distributions of these backgrounds are shown in Fig. [6j These histograms are used to 
model the background contributions in the Dalitz plot fit. 

Using the results of the moments analysis of Sec. [5] as a guide, the nominal Dalitz 
plot fit model for B~ —y D + K~ti~ decays is determined by considering several resonant, 
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Table 5: Signal contributions to the fit model, where parameters and uncertainties are taken 
from Ref. j9j. States labelled with subscript v are virtual contributions. 


Resonance 

Spin 

DP axis 

Model 

Parameters 

£>o(2400)° 

0 

7J1 2 (Dtt) 

RBW 

m = 2318 ± 29 MeV, T = 267 ± 40 MeV 

D*(2460)° 

2 

m 2 (DTT) 

RBW 

Determined from data (see Table [6j) 

£>}( 2760)° 

1 

m 2 (D7i) 

RBW 


Nonresonant 

0 

m 2 (D7i) 

EFF 

Determined from data (see text) 

Nonresonant 

1 

m 2 (Dn) 

FIT' 

£>*(2007)° 

1 

m 2 (D7i) 

RBW 

m = 2006.98 ± 0.15 MeV, V = 2.1 MeV 

Bf 

1 

m 2 (DK) 

RBW 

m = 5325.2 ± 0.4 MeV, T = 0.0 MeV 


nonresonant and virtual amplitudes. Those that do not contribute significantly and that 
do not aid the stability of the fit- are removed. Only natural spin-parity intermediate 
states are considered, as unnatural spin-parity states do not decay to two pseudoscalars. 
The resulting signal model, referred to below as the nominal DP model, consists of the 
seven amplitudes shown in Table [5j three resonances, two virtual resonances and two 
nonresonant terms. Parts of the model are known to be approximations. In particular both 
S- and P-waves in the Dir system are modelled with overlapping broad structures. The 
nominal model gives a better description of the data than any of the alternative models 
considered; alternative models are used to assign systematic uncertainties as discussed in 
Sec. [8l 

The free parameters in the fit are the Cj terms introduced in Eq. (|8]), with the real 
and imaginary parts of these complex coefficients determined for each amplitude in the 
fit model. The £>2(2460)° component-, as the reference amplitude, is the exception with 
real and imaginary parts fixed to 1 and 0, respectively. Fit- fractions and interference fit 
fractions are derived from these free parameters, as are the magnitudes and phases of the 
complex coefficients. Statistical uncertainties for the derived parameters are calculated 
using large samples of simulated pseudoexperiments to ensure that non-trivial correlations 
are accounted for. Several other parameters are also determined from the fit as described 
below. 

In Dalitz plot fits it is common for the minimisation procedure to find local minima of 
the likelihood function. To find the global minimum, the fit is performed many times using 
randomised starting values for the complex coefficients. In addition to the global minimum 
of the likelihood, corresponding to the results reported below, several additional minima 
are found. Two of these have negative log-likelihood (NLL) values close to that of global 
minimum. The main differences between secondary minima and the global minimum are 
the interference patterns in the Dtx S- and P-waves, as shown in App. [Aj 

The shape parameters, defined in Eq. (15), for the nonresonant components are 
determined from the fit- to data to be 0.36 ± 0.03 GeV -2 and 0.36 ± 0.04 GeV -2 for the 
S-wave and P-wave, respectively, where the uncertainties are statistical only. The mass and 
width of the £>2(2460)° resonance are determined from the fit t-o improve the fit quality. 
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Table 6: Masses and widths determined in the fit to data, with statistical uncertainties only. 


Resonance 

Mass (MeV) 

Width (MeV) 

D*(2460)° 

D}(2760)° 

2464.0 ± 1.4 
2781 ± 18 

43.8 ±2.9 

177 ± 32 


Table 7: Complex coefficients and fit fractions determined from the Dalitz plot fit. Uncertainties 
are statistical only. 


Resonance 

Fit fraction (%) 

Real part 

Isobar model coefficients 
Imaginary part Magnitude 

Phase 

£>5(2400)° 

8.3 ±2.6 

-0.04 ±0.07 

-0.51 ±0.07 

0.51 ±0.09 

-1.65 ±0.16 

D* (2460)° 

31.8 ± 1.5 

1.00 

0.00 

1.00 

0.00 

DR2760) 0 

4.9 ± 1.2 

-0.32 ±0.06 

-0.23 ±0.07 

0.39 ±0.05 

-2.53 ±0.24 

S-wave nonresonant 

38.0 ±7.4 

0.93 ±0.09 

-0.58 ±0.08 

1.09 ±0.09 

-0.56 ±0.09 

P-wave nonresonant 

23.8 ±5.6 

-0.43 ±0.09 

0.75 ±0.09 

0.87 ±0.09 

2.09 ±0.15 

D*(2007)° 

7.6 ±2.3 

0.16 ±0.08 

0.46 ±0.09 

0.49 ±0.07 

1.24 ±0.17 

b: 

3.6 ± 1.9 

-0.07 ±0.08 

0.33 ±0.07 

0.34 ±0.06 

1.78 ±0.23 

Total fit fraction 

118.1 






Since the mass and width of the D}(2760)° state have not been precisely determined by 
previous experiments, these parameters are also allowed to vary in the fit. The masses 
and widths of the 0)° and D}(2760)° are reported in Table [6} 

The spin of the _D}(2760)° state has not been determined previously. Fits are performed 
with all values up to 3, and spin 1 is found to be preferred with changes relative to the 
spin 0, 2 and 3 hypotheses of 2ANLL = 37.3,49.5 and 48.2 units, respectively. For 
comparison, the value of 2ANLL obtained from a fit with the 2760)° state excluded 
is 75.0 units. The alternative models discussed in Sec. [8] give very similar values and 
therefore do not affect the conclusion that the D}(2760)° state has spin 1. 

The values of the complex coefficients and fit fractions returned by the fit are shown in 
Table [7J Results for the interference fit fractions are given in App. |Bj The total fit fraction 
exceeds unity mostly due to interference between the Dq( 2400) 0 and S-wave nonresonant 
contributions. 

The consistency of the fit model and the data is evaluated in several ways. Numerous 
one-dimensional projections (including several shown below and those shown in Sec. [5]) 
show good agreement. A two-dimensional y 2 value is determined by comparing the data 
and the fit model in 100 equally populated bins across the SDP. The pull, i.e. the difference 
between the data and fit model divided by the uncertainty, is shown with this SDP 
binning in Fig. [7j The y 2 value obtained is found to be within the bulk of the distribution 
expected from simulated pseudoexperiments. Other unbinned fit quality tests 48 also 
show acceptable agreement between the data and the fit model. 

Figure [8] shows projections of the nominal fit model and the data onto m(Dn), m(DK) 
and m(Kn). Zooms are provided around the resonant structures on iti(Dtt) in Fig. [9] 
Projections of the cosine of the helicity angle of the Dn system are shown in Fig. 10 Good 
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Figure 7: Differences between the data SDP distribution and the fit model across the SDP, in 
terms of the per-bin pull. 


agreement is seen between the data and the fit model. 


8 Systematic uncertainties 

Sources of systematic uncertainty are divided into two categories: experimental and 
model uncertainties. The sources of experimental systematic uncertainty are: the signal 
and background yields in the signal region; the SDP distributions of the background 
components; the efficiency variation across the SDP; possible fit bias. The considered 
model uncertainties are: the fixed parameters in the amplitude model; the addition or 
removal of marginal amplitudes; the choice of models for the nonresonant contributions. 
The systematic uncertainties from each source are combined in quadrature. 

The signal and background yields in the signal region are determined from the fit to 
the B candidate invariant mass distribution, as described in Sec. 4.1. The uncertainty on 


each yield (including systematic uncertainty evaluated as in Sec. 4.3) is calculated, and 
the yields varied accordingly in the DP fit. The deviations from the nominal DP fit result 
are assigned as systematic uncertainties. 

The effect of imperfect knowledge of the background distributions over the SDP 
is tested by varying the histograms used to model the shapes within their statistical 
uncertainties. For ZT*) + 7r _ 7r _ decays the ratio of the D* + and D + contributions is varied. 
Where applicable, the reweighting of the SDP distribution of the simulated samples is 
removed. 

The uncertainty related to the knowledge of the variation of efficiency across the 
SDP is determined by varying the efficiency histograms before the spline fit is performed. 
The central bin in each cell of 3 x 3 bins is varied by its statistical uncertainty and the 
surrounding bins in the cell are varied by interpolation. This procedure accounts for 
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Figure 8: Projections of the data and amplitude fit onto (a) m(Dir), (c) m(DK ) and (e) m(K tt), 
with the same projections shown in (b), (d) and (f) with a logarithmic y -axis scale. Components 
are described in the legend. 
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Figure 9: Projections of the data and amplitude fit onto m(Dir) in (a) the threshold region, 
(b) the Z?2(2460)° region and (c) the 2760)° region. Components are as shown in Fig. |8j 


possible correlations between the bins, since a systematic effect on a given bin is likely also 
to affect neighbouring bins. The effects on the DP fit results are assigned as systematic 
uncertainties. An additional systematic uncertainty is assigned by varying the binning 
scheme of the control sample used to determine the PID efficiencies. 

Systematic uncertainties related to possible intrinsic fit bias are investigated using 
an ensemble of pseudoexperiments. Differences between the input and fitted values from 
the ensemble for the fit- parameters are found to be small. Systematic uncertainties are 
assigned as the sum in quadrature of the difference between the input and output values 
and the uncertainty on the mean of the output value determined from a fit to the ensemble. 

Systematic uncertainties due to fixed parameters in the fit model are determined 
by varying the parameters within their uncertainties and repeating the fit. The fixed 
parameters considered are the mass and width of the Dg(2400)° resonance and the 
Blatt-Weisskopf barrier radius, Tbw- The mass and width are varied by the uncertainties 
shown in Table [ 5 ] and the barrier radius is varied between 3 and 5 GeV -1 jHj. For each fit 
parameter, the difference compared to the nominal fit model is assigned as a systematic 
uncertainty for each source. 

The marginal B*° component is removed from the model and the changes in the 
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Figure 10: Projections of the data and amplitude fit onto the cosine of the helicity angle for the 
Dtt system in (a) the threshold region, (b) the Z?2(2460) 0 region and (c) the D\{ 2760)° region. 
Components are as shown in Fig. [8l 


other parameters are assigned as the systematic uncertainties. Dalitz plot analysis of 
B® —y D°K + 7i~ revealed that a structure at m(D°K + ) ~ 2.86 GeV has both spin 1 and 
spin 3 components j4||5|. Although there is no evidence for a spin 3 resonance in this 
analysis, the excess at m{D + n~) ~ 2.76 GeV could have a similar composition. A putative 
1^3(2760) resonance is added to the fit model, and the effect on the other parameters is 
used to assign systematic uncertainties. 

The EFF lineshapes used to model the nonresonant S- and P-wave contributions are 
replaced by a power-law model and the change in the fit parameters used as a systematic 
uncertainty. The dependence of the results on the effective pole mass description of Eq. (14) 
that is used for the virtual resonance contributions is found by using a fixed width in 


eff 


Eq. (12), removing the dependency on m\ 

The total experimental and model systematic uncertainties for fit fractions and complex 
coefficients are summarised in Tables [8] and [9j respectively. The contributions for the 
fit fractions, masses and widths are broken down in Tables 10 and 11 Similar tables 


summarising the systematic uncertainties on the interference fit fractions are given in 
App-0 The largest source of experimental systematic uncertainty on the fit fractions is 


22 
























































































Table 8: Experimental systematic uncertainties on the fit fractions and complex amplitudes. 





Isobar model coefficients 


Resonance 

Fit fraction (%) 

Real part 

Imaginary part 

Magnitude 

Phase 

Do(2400)° 

0.6 

0.03 

0.02 

0.02 

0.06 

D*(2460)° 

0.9 

- 

- 

- 

- 

£*(2760)° 

0.4 

0.03 

0.03 

0.01 

0.08 

S-wave nonresonant 

1.5 

0.03 

0.03 

0.02 

0.04 

P-wave nonresonant 

2.1 

0.03 

0.05 

0.03 

0.05 

£*(2007)° 

1.3 

0.03 

0.04 

0.04 

0.07 

B* v 

0.9 

0.22 

0.02 

0.03 

0.11 

Table 9: Model uncertainties on 

the fit fractions and complex amplitudes. 





Isobar model coefficients 


Resonance 

Fit fraction (%) 

Real part 

Imaginary part 

Magnitude 

Phase 

£*(2400)° 

1.9 

0.28 

0.13 

0.15 

0.51 

£^(2460)° 

1.4 

- 

- 

- 

- 

£>*(2760)° 

0.9 

0.03 

0.03 

0.03 

0.08 

S-wave nonresonant 

10.8 

0.17 

0.15 

0.20 

0.11 

P-wave nonresonant 

3.7 

0.34 

0.68 

0.12 

0.95 

£*(2007)° 

1.5 

0.56 

0.77 

0.05 

0.60 

B* v 

1.6 

0.09 

0.08 

0.07 

0.27 


due to the efficiency variation. For the model uncertainty on the fit fractions, the addition 
and removal of marginal components and variation of fixed parameters dominate. In 
general, the model uncertainties are larger than the experimental systematic uncertainties 
for the fit fractions and the masses and widths. 

Several cross-checks are performed to confirm the stability of the results. The data 
sample is divided into two parts depending on the charge of the B candidate, the polarity 
of the magnet and the year of data taking. Selection effects are also checked by varying 
the requirement on the neural network output variable and the PID criteria applied to the 
bachelor kaon. A fit is performed for each of the subsamples individually and each is seen 
to be consistent with the default fit results, although in some cases one of the secondary 
minima described in App. [A] becomes the preferred solution. To cross-check the amplitude 
model, the fit is repeated many times with an extra resonance with fixed mass, width and 
spin included in the model. All possible mass and width values, and spin up to 3, were 
considered. None of the additional resonances are found to contribute significantly. 
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Table 10: Breakdown of experimental systematic uncertainties on the fit fractions (%) and masses 
(MeV) and widths (MeV). 



Nominal 

S/B frac. Eff. 

Bkg. 

Fit bias 

Total 

£q(2400)° 

8.3 ±2.6 

0.2 

0.5 

0.1 

0.3 

0.6 

H 2 (2460)° 

31.8 ±1.5 

0.2 

0.8 

0.0 

0.2 

0.9 

£^(2760)° 

4.9 ±1.2 

0.2 

0.2 

0.1 

0.2 

0.3 

S-wave nonresonant 

38.0 ±7.4 

0.7 

0.5 

0.4 

1.2 

1.5 

P-wave nonresonant 

23.8 ±5.6 

1.0 

1.6 

0.7 

0.5 

2.1 

£);(2007) 0 

7.6 ±2.3 

0.7 

1.0 

0.3 

0.3 

1.3 

K 

3.6 ±1.9 

0.3 

0.3 

0.2 

0.8 

0.9 

m{D*( 2460)°) 

2464.0 ± 1.4 

0.1 

0.1 

0.0 

0.2 

0.2 

T{D* (2460)°) 

43.8 ±2.9 

0.3 

0.3 

0.0 

0.4 

0.6 

m(£*(2760)°) 

2781 ± 18 

1 

4 

0 

2 

6 

T (7^(2760)°) 

177 ± 32 

3 

1 

2 

5 

7 

Table 11: Breakdown of model uncertainties on the fit fractions (%) and masses (MeV) and 

widths (MeV). 








Nominal Add / rem 

Alt. models 

Fixed params Total 

£>*(2400)° 

8.3 ±2.6 

2.0 

0.1 


0.2 

2.0 

.Do (2460)° 

31.8 ±1.5 

1.3 

0.2 


0.4 

1.4 

£*(2760)° 

4.9 ± 1.2 

0.8 

0.1 


0.3 

0.9 

S-wave nonresonant 

38.0 ± 7.4 

4.8 

4.5 


5.4 

10.8 

P-wave nonresonant 

23.8 ±5.6 

2.6 

2.1 


3.0 

3.7 

£*(2007)° 

7.6 ±2.3 

0.6 

0.1 


1.4 

1.5 

B* v 

3.6 ±1.9 

0.7 

1.0 


1.1 

1.6 

m(D*( 2460)°) 2464.0 ± 1.4 

0.5 

0.1 


0.1 

0.5 

r (£3(2460)°) 

43.8 ±2.9 

0.8 

1.4 


0.6 

1.7 

m(£*(2760)°) 

2781 ± 18 

6 

6 


1 

11 

r(£^( 2760 )°) 

177 ± 32 

16 

9 


1 
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9 Results and summary 


The results for the complex coefficients are reported in Tables [12] and [TTT| in terms of real and 
imaginary parts and of magnitudes and phases, respectively. The results for the fit fractions 
are given in Table p~4| and the results for the interference fit fractions are given in App.[B} The 
fit fractions for resonant contributions are converted into quasi-two-body product branching 
fractions by multiplying by B(B~ —> D + K~tt~) = (7.31 ± 0.19 ± 0.22 ± 0.39) x 10 -5 , as 
determined in Sec. |4.4 These product branching fractions are shown in Table 15 they 
cannot be converted into absolute branching fractions because the branching fractions for 


24 



























Table 12: Results for the complex amplitudes and their uncertainties. The three quoted errors 
are statistical, experimental systematic and model uncertainties, respectively. 


Resonance 

Isobar model coefficients 

Real part Imaginary part 

£>*(2400)° 

±>*(2460)° 

£>*(2760)° 

-0.04 ±0.07 ±0.03 ±0.28 

1.00 

-0.32 ±0.06 ±0.03 ±0.03 

-0.51 ±0.07 ±0.02 ±0.13 

0.00 

-0.23 ±0.07 ±0.03 ±0.03 

S-wave nonresonant 
P-wave nonresonant 

0.93 ±0.09 ±0.03 ±0.17 
-0.43 ±0.09 ±0.03 ±0.34 

-0.58 ±0.08 ±0.03 ±0.15 

0.75 ±0.09 ±0.05 ±0.68 

£>*(2007)° 

B* v 

0.16 ±0.08 ±0.03 ±0.56 

-0.07 ±0.08 ±0.22 ±0.09 

0.46 ±0.09 ±0.04 ±0.77 
0.33 ±0.07 ±0.02 ±0.08 


Table 13: Results for the complex amplitudes and their uncertainties. The three quoted errors 
are statistical, experimental systematic and model uncertainties, respectively. 


Resonance 

Isobar model coefficients 

Magnitude Phase 

-Dg(2400)° 

±>*(2460)° 

±>*(2760)° 

0.51 ±0.09 ±0.02 ±0.15 

1.00 

0.39 ±0.05 ±0.01 ±0.03 

-1.65 ±0.16 ±0.06 ±0.50 

0.00 

-2.53 ±0.24 ±0.08 ±0.08 

S-wave nonresonant 
P-wave nonresonant 

1.09 ±0.09 ±0.02 ±0.20 

0.87 ±0.09 ±0.03 ±0.11 

-0.56 ±0.09 ±0.04 ±0.11 

2.09 ±0.15 ±0.05 ±0.95 

±>*(2007)° 

b: 

0.49 ±0.07 ±0.04 ±0.05 
0.34 ±0.06 ±0.03 ±0.07 

1.24 ±0.17 ±0.07 ±0.60 
1.78 ±0.23 ±0.11 ±0.27 


Table 14: Results for the fit fractions and their uncertainties (%). The three quoted errors are 
statistical, experimental systematic and model uncertainties, respectively. 


Resonance 

Fit fraction 


±>*(2400)° 

8.3 ± 2.6 ± 0.6 ± 

1.9 

±>*(2460)° 

31.8 ± 1.5 ± 0.9 ± 

1.4 

±>*(2760)° 

4.9 ± 1.2 ± 0.3 ± 

0.9 

S-wave nonresonant 

38.0 ±7.4 ±1.5 ±10.8 

P-wave nonresonant 

23.8 ± 5.6 ± 2.1 ± 

3.7 

H*(2007)° 

7.6 ± 2.3 ± 1.3 ± 

1.5 

b : 

3.6 ± 1.9 ± 0.9 ± 

1.6 


the resonance decays to D + n are unknown. 


25 



Table 15: Results for the product branching fractions B(B~ — > RK ~) x B(R —> D + ir ~) (10~ 6 ). 
The four quoted errors are statistical, experimental systematic, model and inclusive branching 
fraction uncertainties, respectively. 


Resonance 

Branching fraction 

£>o(2400)° 

6.1 ±1.9 ±0.5 ±1.4 ±0.4 

£>*(2460)° 

23.2 ±1.1 ±0.6 ±1.0 ±1.6 

£>*(2760)° 

3.6 ±0.9 ±0.3 ±0.7 ±0.2 

S-wave nonresonant 

27.8 ±5.4 ±1.1 ±7.9 ±1.9 

P-wave nonresonant 

17.4 ±4.1 ± 1.5 ±2.7 ± 1.2 

£*(2007)° 

5.6 ±1.7 ±1.0 ±1.1 ±0.4 

B* 

2.6 ±1.4 ±0.6 ±1.2 ±0.2 


The masses and widths of the T>2(2460)° and D\{ 2760)° are determined to be 


m(£>;(2460)°) 

T(£>;(2460)°) 

m(£*(2760)°) 

r(£>*(2760)°) 


(2464.0 ± 1.4 ± 0.5 ± 0.2) MeV , 
(43.8 ± 2.9 ± 1.7 ±0.6) MeV, 
(2781 ± 18 ± 11 ± 6) MeV, 
(177 ± 32 ± 20 ± 7) MeV, 


where the three quoted errors are statistical, experimental systematic and model uncer¬ 
tainties, respectively. The results for the £>2(2460)° are within 2 er of the world average 
values |9|. The mass of the £>*(2760)° resonance is similarly consistent with previous 
measurements. The measured width of this state is larger than previous measurements by 
2 to 3 times the uncertainties. Future studies based on much larger data samples will be 
required to better understand these states. 

The measurement of B(B~ —> D + K~n~) corresponds to the first observation of this 
decay mode. Therefore, the resonant contributions to the decay are also first observations. 
The significance of the B~ —> £>^(2760)°/! ~ observation is investigated by removing the 
corresponding resonance from the DP model. A fit without the £>^(2760)° component 
increases the value of 2ANLL by 75.0 units, corresponding to a high statistical significance. 
Only the systematic effects due to uncertainties in the DP model could in principle 
significantly change the conclusion regarding the need for this resonance. However, in 
alternative DP models where a Dn resonance with spin 3 is added and where the B* 
contribution is removed, the shift in 2ANLL remains above 50 units. The alternative models 
also do not significantly impact the level at which the 2760)° state is preferred to be 
spin 1. Therefore, these results represent the first observation of the B~ —> £]“(2760) 0 /l _ 
and the measurement of the spin of the £^(2760)° resonance. 

In summary, the B~ —>■ D + K~tt~ decay has been observed in a data sample corre¬ 
sponding to 3.0 ftV 1 of pp collision data recorded by the LHCb experiment. An amplitude 
analysis of its Dalitz plot distribution has been performed, in which a model containing 
resonant contributions from the £>q( 2400)°, £>2(2460)° and D\{ 2760)° states in addition 
to both S-wave and P-wave nonresonant amplitudes and components due to virtual 
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D*{ 2007)° and B*° resonances was found to give a good description of the data. The 
B~ —y £)2(2460)°/i _ decay may in future be used to determine the angle 7 of the CKM 
unitarity triangle. The results provide insight into the spectroscopy of charm mesons, and 
demonstrate that further progress may be obtained with Dalitz plot analyses of larger 
data samples. 
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A Secondary minima 


The results, in terms of fit- fractions and complex coefficients, corresponding to the two 
secondary minima discussed in Sec. [7] are compared to those of the global minimum 
in Table [16} The main difference between the global and secondary minima is in the 


interference pattern in the Dtt P-waves, while the third minimum exhibits a different- 
interference pattern in the Dtt S-wave than the global minimum and has a very large total 
fit fraction due to strong destructive interference. 


Table 16: Results for the fit fractions and complex coefficients for the secondary minima with 
2NLL values 2.8 and 3.3 units greater than that of the global minimum of the NLL function. 


Resonance 

2ANLL 

Fit fraction (%) 

0 2.8 3.3 

0 

Real part 
2.8 

3.3 

Imaginary part 

0 2.8 3.3 

Magnitude 

0 2.8 3.3 

0 

Phase 

2.8 

3.3 

n;(2400) u 

8.3 

9.6 

84.4 

-0.04 

-0.03 

-1.38 

-0.51 

-0.55 

-0.72 

0.51 

0.55 

1.56 

-1.65 

-1.62 

-2.66 

£>J;(2460)° 

31.8 

31.5 

34.9 


1.00 



0.00 



1.00 



0.00 


£>1(2760)° 

4.9 

4.6 

5.5 

-0.32 

-0.30 

-0.30 

-0.23 

-0.24 

-0.26 

0.39 

0.38 

0.40 

-2.53 

-2.46 

-2.42 

S-wave nonresonant 

38.0 

36.2 

4.6 

0.93 

0.89 

-0.33 

-0.58 

-0.60 

0.15 

1.09 

1.07 

-0.36 

-0.56 

-0.59 

2.71 

P-wave nonresonant 

23.8 

22.6 

-31.9 

-0.43 

0.83 

-0.84 

0.75 

0.15 

0.45 

0.87 

0.85 

0.96 

2.09 

2.96 

2.65 

£>;(2007)“ 

7.6 

7.1 

11.9 

0.16 

-0.38 

-0.28 

0.46 

-0.29 

-0.51 

0.49 

0.48 

0.58 

1.24 

-2.49 

-2.07 

B* 

3.6 

1.0 

25.0 

-0.07 

-0.16 

-0.31 

0.33 

0.09 

0.79 

0.34 

0.18 

0.85 

1.78 

2.61 

1.94 

Total fit fraction 

118.1 

112.6 

198.3 



B Results for interference fit fractions 


The central values and statistical errors for the interference fit- fractions are shown in 


Table 17 The experimental systematic and model uncertainties are given in Tables 18 


and 19 The interference fit fractions are common to both trigger subsamples. 


Table 17: Interference fit fractions (%) and statistical uncertainties. The amplitudes are: (To) 
£>*(2007)°, (Ti) £>*(2400)°, (T 2 ) £>^(2460)°, (T 3 ) £>*(2760)°, (T 4 ) B*, (T 5 ) nonresonant S-wave, 
(A§) nonresonant P-wave. The diagonal elements are the same as the conventional fit fractions. 



Ao 

Ti 

A-2 

T3 

A 4 

T5 

Te 

To 

7.6 ± 2.3 

0.0 ±0.0 

0.0 ± 0.0 

2.4 ±0.9 

4.8 ± 1.3 

0.0 ±0.0 

-14.2 ±5.3 

Ti 


8.3 ±2.6 

0.0 ± 0.0 

0.0 ±0.0 

— 1.6 ± 0.7 

18.1 ±2.6 

0.0 ±0.0 

^2 



31.8 ±1.5 

0.0 ±0.0 

-2.3 ±0.6 

0.0 ±0.0 

0.0 ±0.0 

T3 




4.9 ±1.2 

2.0 ±0.8 

0.0 ±0.0 

-9.6 ±2.9 

a 4 





3.6 ± 1.9 

-6.7 ±2.3 

-11.1 ± 3.6 

T5 






38.0 ± 7.4 

0.0 ±0.0 

To 







23.8 ±5.6 
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Table 18: Experimental systematic uncertainies on the interference fit fractions (%). The 
amplitudes are: (M 0 ) £>*( 2007 )°, (M 4 ) £>5(2400)°, (M 2 ) £> 2 (2460)°, (M 3 ) D\{ 2760)°, (M 4 ) B*, 
(M 5 ) nonresonant S-wave, (A§) nonresonant P-wave. The diagonal elements are the same as the 
conventional fit fractions. 



Aq 


A 2 

A 3 

7U 

A 5 

Ae 

Aq 

1.3 

0.0 

0.0 

0.4 

0.6 

0.0 

2.6 

A\ 


0.6 

0.0 

0.0 

0.4 

0.6 

0.0 

A 2 



0.9 

0.0 

0.3 

0.0 

0.0 

A 3 




0.4 

0.2 

0.0 

0.7 

a 4 





0.9 

1.1 

1.2 

A 5 






1.5 

0.0 

Aq 







2.1 


Table 19: Model systematic uncertainies on the interference fit fractions (%). The amplitudes are: 
(M 0 ) £>*(2007)°, (Mi) £>q(2400)°, (M 2 ) £>|(2460)°, (M 3 ) £>^(2760)°, (M 4 ) B*, (M 5 ) nonresonant 
S-wave, (M 6 ) nonresonant P-wave. The diagonal elements are the same as the conventional fit 
fractions. 



Aq 

A 1 

A 2 

^3 

A 4 

A 5 

Aq 

Ao 

1.5 

0.0 

0.0 

0.1 

1.4 

0.0 

1.1 

A 1 


1.9 

0.0 

0.0 

1.7 

4.8 

0.0 

A 2 



1.4 

0.0 

0.5 

0.0 

0.0 

A 3 




0.9 

0.6 

0.0 

3.4 

a 4 





1.6 

2.8 

0.4 

A 5 






10.8 

0.0 

Aq 







3.67 
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